Noise thresholds for optical quantum computers 
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In this paper we numerically investigate the fault-tolerant threshold for optical cluster-state quan- 
tum computing. We allow both photon loss noise and depolarizing noise (as a general proxy for 
all local noise), and obtain a threshold region of allowed pairs of values for the two types of noise. 
Roughly speaking, our results show that scalable optical quantum computing is possible for photon 
loss probabilities < 3 x 10~ 3 , and for depolarization probabilities < 10 -4 . 
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Optical systems are promising candidates for quantum 
computation, due to their long decoherence times, accu- 
rate single-qubit gates, and relatively efficient readout. 
A scheme for optical quantum computing has been sug- 
gested by Knill, LafLamme and Milburn (KLM) pj, and 
the basic elements of that scheme experimentally demon- 
strated 

HH 110 Unfortunately, KLM requires tens 
of thousands of optical elements to achieve a single en- 
tangling gate operating with high probability. A recent 
proposal jg (c.f. 0) combines elements of KLM with the 
cluster-state model of quantum computation to re- 
duce the complexity by many orders of magnitude. This 
scheme has been simplified || to require only tens of op- 
tical elements per logical gate. A recent experiment 
demonstrated simple optical cluster state computing. 

Our paper investigates the effect of noise on the optical 
cluster state proposals. In the standard quantum circuit 
model, the noise threshold theorem shows that provided 
the amount of noise per elementary operation is below 
the threshold, scalable quantum computation is possible. 
Unfortunately, this result does not apply directly to opti- 
cal cluster states, as the cluster state model is fundamen- 
tally different from the circuit model. However, [Til [T^ 
(c.f. E3) established the existence of a threshold for clus- 
ters, without obtaining a value, while 

El 

argued that 

in a certain noise model the cluster threshold is no more 
than an order of magnitude worse than the circuit thresh- 
old. This latter work is not directly relevant to optical 
clusters, since it uses deterministic entangling gates, and 
does not include any process analogous to photon loss. 

We use numerical simulations to find the noise thresh- 
old. Our analysis is tailored to the dominant sources of 
noise in optical quantum computing, including the non- 
determinism of the optical entangling gates, photon loss, 
and depolarizing noise. We therefore obtain a thresh- 
old region of noise parameters for which scalable optical 
quantum computing is possible. Our protocol is complex, 
and we omit some details; full details will appear in [l5j . 

A prior work 16] has calculated a threshold for optical 
quantum computation when the only source of noise is 
photon loss. In real experiments other noise sources such 
as dephasing are also present, and protocols such as 
will actually amplify the effects of such noise at the en- 



coded level. By contrast, our protocol protects against 
both photon loss and depolarizing noise, and by stan- 
dard fault-tolerance results thus automatically protects 
against arbitrary local noise, including dephasing (in any 
basis), amplitude damping, etc. 

Introductions to cluster-state computation may be 
found in 0, and we assume familiarity with the 
model. An important element in the model are the Pauli 
"byproduct" operators, known functions of the measure- 
ment results, which are used to correct the state when 
the computation concludes. We call the tensor product 
of these byproduct operators the Pauli frame, and it is 
updated after measurements of cluster qubits according 
to a set of propagation rules, described, e.g., in [l7| . 

Our approach to optical cluster-state computation is 
based most closely on [jj. We use the polarization of a 
single photon to encode a single qubit, and build clus- 
ters up using fusion gates ("type I fusion gates" in |9(), 
which, when applied to two cluster qubits either (a) fuse 
the qubits into a single cluster qubit, which occurs with 
probability |; or (b) measure both qubits in the compu- 
tational basis, also with probability \. 

Available physical resources and noise: Our re- 
sources are: (1) a source of polarization entangled Bell 
pairs; (2) single-qubit gates, effected using linear optics, 
as in KLM; (3) efficient polarization-discriminating pho- 
ton counters capable of distinguishing 0, 1 and 2 photons; 
(4) fusion gates, built from beamsplitters and photon 
counters; and (5) quantum memory gates. We assume 
all these elements take the same amount of time, and 
describe our circuit as a sequence of such time steps. 

Our noise model includes a parameter 7 representing 
the probability per qubit per time step of photon loss. 
We assume this probability is independent of the state 
of the qubit, and that photon loss occurs after Bell-state 
preparation, and before memory, single-qubit and fusion 
gates; for two-qubit operations we assume photon loss 
occurs independently for both qubits. Our noise model 
also includes a depolarizing parameter e. Depolarizing 
noise affects physical operations as follows: (1) after Bell- 
state preparation or before a fusion gate the two qubits 
are collectively depolarized, i.e., with probability 1 — e 
nothing happens, while with probabilities e/15 we apply 
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each of the 15 non-identity Pauli operators IX, XX etc; 
and (2) before memory, single-qubit and measurement 
gates the qubit is depolarized with parameter e. 

Additional noise sources that may effect real imple- 
mentations include dark counts and dcphasing. However, 
the fault-tolerant protocol we implement automatically 
protects against such noise sources, and we believe the 
threshold results will not qualitatively change. 

Method of simulation: We use the stabilizer for- 
malism to simulate Clifford group operations, which are 
sufficient to simulate error-correction and depolarization. 
In our simulations, rather than working with the state 
directly, we merely keep track of the errors in the state 
when compared with an ideal reference state. We keep 
track of two types of errors: the physical error in the state 
of the cluster, which is represented as a tensor product 
of Pauli operators, and errors in the Pauli frame due to 
erroneous measurement results, which are again a tensor 
product of Pauli operators. Note that physical errors may 
propagate to become Pauli frame errors when qubits with 
physical noise are measured, giving rise to incorrect mea- 
surement results. The rules for propagating both types 
of errors may be computed following, e.g., [17|: see |15j . 

These methods suffice to describe error-correction and 
Pauli-type noise, but not photon loss and fusion gate 
failures. We can use postselection and repetition to ef- 
fectively eliminate photon loss and fusion gate failures, 
whenever those failures do not directly affect the encoded 
data. However, when they do affect the data, another 
approach must be taken. Suppose when fusion gate fail- 
ure occurs the experimenter: (1) Randomizes the local 
Pauli frame of the data qubit; (2) notes the location at 
which the failure occurred, for use in decoding; and (3) 
carries out the rules for propagating the Pauli frame, as 
though the fusion gate had succeeded, and the bond was 
created. The rules for propagating Pauli frame errors 
can be used to show that once the experimenter has ran- 
domized the Pauli frame, it does not make any physical 
difference whether the fusion gate failed or not, and so 
we can treat it as though it succeeded. The remaining er- 
rors are Pauli-type errors, and so can be simulated in the 
standard way. The details (and a discussion of photon 
loss, which is dealt with similarly) appear in [15j. 

Broad picture of fault-tolerant protocol: The 
protocol is split into two parts: (1) a cluster-based simu- 
lation of a variant of Steane's protocol 0|; (2) a deter- 
ministic gate-based protocol, again based on 191. The 
cluster threshold is obtained by concatenating the results 
from a single level of the cluster protocol with multiple 
levels of the deterministic protocol. The idea is to take 
a quantum circuit, build up a fault-tolerant simulation 
through multiple levels of concatenation in the circuit 
model, and then replace the bottom level by a cluster- 
ized simulation of a noisy deterministic gate. 

Microclusters and parallel fusion: Our protocol is 
based on microclusters, star-shaped clusters with a cen- 
tral root node, and attached leaf nodes. Such microclus- 
ters can be created using repeated fusion of Bell pairs. By 



attempting preparation of a large number of microclus- 
ters in parallel and postselecting on successful attempts 
we can build a fc-leaf microcluster in 0(log(fc)) timesteps, 
and consuming 0(k 2 ) Bell pairs, with probability of suc- 
cess arbitrarily close to one. 

Microclusters can be used to ensure that larger clus- 
ters always have multiple leaf nodes. This can be used 
to enhance the probability of fusing two clusters, by at- 
tempting simultaneous fusion gates between adjacent leaf 
nodes of the two clusters. With a probability that goes 
rapidly to one as the number of leaves increases, at least 
one of these fusion gates succeeds, fusing the two clusters 
together. We call this process of using leaves to fuse the 
two clusters with high probability parallel fusion. 

State at the start of a round: The ideal noise-free 
state at the start of any round is the encoded state of 
the data (e.g., in the 7-qubit code), but with a number 
of leaves attached to each code qubit, which are used 
later for parallel fusion. 

At the beginning of the entire trial, we assume the 
input is a noiseless state of this form. We justify this 
assumption on the grounds that the initial state does 
not actually matter, since our goal is to estimate the 
rate per round at which crashes occur in the encoded 
data. Following i we perform warm-up rounds of error 
correction before gathering data on this crash rate, so as 
to avoid transient effects due to the choice of initial state. 

Ancilla creation: Each round of error correction 
involves the creation of multiple verified ancilla states, 
which are used to extract syndrome bits. We illustrate 
this for the 7-qubit code, but the procedure generalizes 
to other codes. We create the ancilla using the cluster: 




(1) 

This is a clusterization of the ancilla creation circuit 
in [19j , with the final column of the cluster correspond- 
ing to an encoded |+). We abridge our notation so that 
touching circles represent connected cluster qubits. The 
cluster is created by first creating an array of microclus- 
ters; the large circles represent root nodes, while the 
smaller circles represent leaves; note that many of the 
leaves are consumed during preparation by fusion and 
parallel fusion, and are not shown. We then use fusion 
and parallel fusion to create the bonds; details appear 
in 15]. We conclude by measuring all qubits in the X 
basis, except the leaves in column 8. To verify the ancilla, 
we postselect on the measurement results of the termi- 
nating qubits in rows 0, 1, 2, 3 all being 0. The resulting 
state is an encoded |+), with each qubit having a number 
of leaves attached for the purpose of parallel fusion. 

The telecorrector: To extract error syndromes we 
interact the data and ancillas using a special cluster state 
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called a telecorrector. Telecorrector-based syndrome ex- 
traction is a variant of Steane's approach (c.f. also the 
related protocol in 20]), and can be thought of as a clus- 
terized version of the circuit: 
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where operations are being performed on encoded qubits, 
0} is fault-tolerant ancilla creation, and the measurement 
is a transversal A basis measurement. 

Telecorrector creation begins with the creation of mul- 
tiple copies of the following state, one copy for each qubit 
in the code being used: 





(3) 

This state can be created in the obvious way using post- 
selected microcluster fusion. The leaves on the left-hand 
end will eventually be used to attach to a single qubit of 
the encoded data using parallel fusion. The leaves and 
root node on the right-hand end will contain the output 
of this round of error-correction, and become the input to 
the next round of error-correction. The remaining leaves 
will be used to fuse to ancilla states. 

Simultaneous with the creation of the state © , we cre- 
ate four verified ancilla states, and fuse the ancillas with 
the leaves on the state © to create the state (illustrated 
as though for a three-qubit code): 
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Next we measure all the shaded qubits in the A basis, 
leaving only the leftmost and rightmost leaves, for later 
use in attaching the data, and future rounds of error cor- 
rection. Applying the propagation rules for the Pauli 
frame, it can be shown that the measurement outcomes 
from the shaded qubits completely determine whether 
the repeated syndrome measurements will agree or not, 
before the state has interacted with the data! We have 
verified this fact both numerically and analytically. 

We take advantage of this remarkable fact by postse- 
lecting on measurement outcomes that ensure this prea- 
greeing syndrome property. We call the resulting state 
the telecorrector. The preagreeing syndrome property 
enables syndrome extraction to be performed more ef- 
ficiently (thus improving the threshold) than in Steane's 
protocol, which extracts many syndromes to ensure that 
some large subset agree. 



Once prepared, we use parallel fusion to attach the 
telecorrector to the data, and then X basis measurements 
to complete this part of the computation. Standard prop- 
agation rules are used to update the Pauli frame, and to 
determine the syndrome extracted from this procedure. 

Decoding: We use a technique for syndrome decoding 
which takes advantage of the experimenter's knowledge 
of the locations of photon loss and fusion gate failures. In 
particular, we use the fact (see p. 467 of |21|1 that a code 
correcting t unlocated errors is able to correct 2t located 
errors. Our technique is a maximum likelihood procedure 
for decoding arbitrary combinations of located and un- 
located errors. All the codes we use are CSS codes with 
the property that decoding of the X and Z errors can be 
performed separately using an identical procedure. The 
A-decoding procedure (for example) has the following 
inputs: the measured A-error syndrome, obtained from 
the vector of total errors of the ancilla measurement out- 
comes; and a list of locations (qubit indices within the 
code block) at which located errors have occurred during 
the round. The outputs of the decoding routine are: a 
list of locations where X flips should be made in order 
to correct the data; and a flag signalling a located crash. 
The located crash flag is set to "true" when different 
patterns of X errors are found to have equal maximum 
likelihood, but differ from each other by a logical A oper- 
ation. The located crash flag is used to improve decoding 
at the next level of concatenation, by identifying encoded 
blocks which are known to have experienced an error. 

Results of the optical cluster simulation: Our 
aim in simulating cluster-based error correction is to es- 
timate the function which maps the input noise parame- 
ters (e, 7) to the logical error rates, or crash rates, defined 
below. Below we describe simulations which estaimate a 
similar function for a deterministic circuit-based proto- 
col, and then combine the results to give the threshold 
curve for cluster-state optical quantum computing. 

At the end of a round of simulated cluster-based error 
correction, we say that the round has caused a located 
crash whenever either the X or Z decoding steps have 
reported a located crash. We define an unlocated crash 
as follows. We take the pattern of Pauli errors on the 
root nodes of the data, and consider the result of a per- 
fect (noise-free) round of correction. If perfect correction 
would result in a pattern on Pauli errors corresponding 
to a non-identity encoded Pauli operation, then we say 
the data has experienced an unlocated crash. 

We performed simulations based on the Golay 23-qubit 
and Steane 7-qubit codes. For each simulation, we chose 
a number of settings for the noise parameters (e, 7), and 
for each we ran a many-trial monte carlo simulation. 
Each trial consisted of two successive rounds of error cor- 
rection, and the outcome of the trial was determined by 
whether the second of the two rounds caused a crash. 
The purpose of the first "warm-up" round is to reduce 
transient effects due to our choice of (noise free) initial 
conditions. Including more than one warm-up round did 
not make a statistically significant change to the results. 
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FIG. 1: Threshold region for optical clusters using the 23- 
qubit Golay code (solid) and 7-qubit Steane code (dashed). 

We tally the outcomes as follows. For all the trials for 
which the first round does not cause a crash, we count: 
(1) the number Njj for which the second round causes an 
unlocated crash but not a located crash, (2) the number 
Nl for which the second round causes a located crash, 
and (3) the number Nn for which no crashes occur. The 
unlocated and located crash rates E and T are estimated 
as E = M N " AT and T = . r . U T L , , T , We calculate 
these crash rates for a variety of input noise parameters, 
and use weighted least-squares fitting to fit polynomials 
E(e,j) and r(e,7) representing the general behaviour of 
the crash rates; these polynomials were in good agree- 
ment with the qualitative theory of fault-tolerance. 

Concatenation: Under k layers of concatenation, our 
protocol is effectively equivalent to doing k — 1 concate- 
nated levels of an ordinary deterministic gate-based fault- 
tolerance protocol, and then replacing the elements at the 
lowest level by cluster-based equivalents with just a sin- 
gle level of encoding. To understand the behaviour of 



the concatenated protocol, we therefore also did simu- 
lations of a deterministic fault-tolerant protocol. These 
simulations followed , but incorporated a determinis- 
tic version of telecorrection, and the separate treatment 
of unlocated and located errors. The appropriate noise 
model has two noise parameters, (p,q), representing, re- 
spectively, the rate of located and unlocated Pauli er- 
rors, corresponding to located and unlocated crashes at 
the next lowest level of concatenation. The results of our 
simulations suggest that these crashes may be accurately 
modelled as independent X and Z errors, with Y errors 
suppressed, and so this is the noise model we adopt. We 
use least squares fitting to estimate polynomials P(p, q) 
and Q(p,q), where P and Q are the rates for located 
and unlocated crashes at a single level of encoding in the 
cluster-based simulations. 

Results and conclusion: Define maps / : (e, 7) — » 
{E,T) and g : (p, q) — > (P,Q). Then the located and 
unlocated crash rates after k levels of concatenation may 
be estimated by computing x ) ° /)(e,7). Provided 
this tends to (0, 0) as k — > 00 we are inside the thresh- 
old region. Fig. 2] illustrates the threshold region, and 
shows that the 23-qubit code gives a considerably bet- 
ter threshold than the 7-qubit code; resource usage will 
be discussed in 0. Both codes give thresholds worse 
than the best known circuit thresholds |2j}, but the re- 
sults are encouraging given the non-deterministic nature 
of the optical entangling operations. 
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